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Abstract 



Positronium decay rates are computed at the one-loop level, using convolution- 
type factorized amplitudes. The dynamics of this factorization is probed 
with dispersion relations, showing that unallowed approximations are usu- 
ally made, and some O (a 2 ) corrections missed. Further, we discuss the 
relevance of the Schrodinger wavefunction as the basis for perturbative 
calculations. Finally, we apply our formalism to the parapositronium two- 
photon decay. 
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1 Introduction 

Positronium is a bound state of electron and positron. In this paper, we will be 
interested in the singlet state, parapositronium, whose decay rate into 27 has 
been precisely measured JTJ: 

T exp (p-Ps -► 2 7 ) = (7.9909 ± 0.0017) x 10 9 sec" 1 (1) 

The corresponding theoretical predictions which include perturbative QED 
corrections to a non-relativistic treatment of the bound state wavefunction have 
been computed also with high accuracy (see for example [2], @, @): 

r {p-Ps -» 27) 



with a the fine structure constant and m the electron mass. As can be ob- 
served, agreement between theory and experiment is good. However, the non- 
logarithmic O (a 2 ) corrections, which have been obtained only recently, are not 
yet accessible experimentally. 

Positronium is a test ground for bound state treatment in Quantum Field 
Theory. The first try dates back to the 40's, with decay rates expressed through 
a factorized formula || 

r (p-Ps - 2 7 ) = |0 O | 2 • (4v rel a (e+e- -> 27))^^ 

with <j) the Schrodinger positronium wavefunction at the origin, a (e + e~ — ► 27) 
the total cross section for e + e~ — > 27 and iv e i the relative velocity of e + and e~ 
in their center of mass frame. Since then, more sophisticated decay amplitudes 
have been constructed, and systematic procedures for calculating corrections 
have been developed. However, the basic factorization of the bound state dy- 
namics from the annihilation process has remained as a basic postulate. For 
low order corrections, this approximation is unquestionable, but for O (a 2 ) cor- 
rections, factorization has to be tested. Indeed, non-perturbative phenomena 
responsible for the off-shellness of the electron and positron inside the positro- 
nium are of O (a 2 ) . In other words, to get a sensible theoretical prediction at 
O (a 2 ), one must carefully analyze how binding energy effects enter the general 
factorization approach. 

In the present paper, we propose a systematic procedure for factorizing the 
bound state dynamics from the annihilation process. From a fully relativistic 
model at the one-loop level, where off-shell constituents appear, we will recover 
the standard factorized amplitude used in the literature. Most importantly, we 
will show that those standard formulas involve some unnecessary approxima- 
tions, and we will remove them. Since our derivation relies on well-established 
techniques of quantum field theory, we conclude that some O (a 2 ) corrections 
have been forgotten. Our derivation of the factorized formula is particularized 
to the parapositronium two-photon decay for the sake of definiteness, but it is 
completely general, so equally valid for orthopositronium decays. 
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a 5 m 



— - + 2a 2 log- + 1.75 30 - - — log 2 - 
4/7T a V 7r / 2-7T a 



7.98950(2) x 10 9 sec _1 



After those general considerations, an alternative and equivalent factorized 
amplitude is found for parapositronium, which allows some further analyses 
of the factorized formula. Using this last form, we compute lowest order plus 
binding energy corrections for parapositronium decay. Those calculations are 
carried with two different input forms for the bound state wavefunction; first 
the well-known Schrodinger momentum wavefunction, and then an improved 
form for this wavefunction. Finally, we comment on those two forms, and 
point towards the fact that it is possibly not the usual simple Schrodinger 
wavefunction, but rather the improved one that should be used, due to the 
dynamics of the factorization. 



2 The Loop Model for Positronium Decay 

In this section, we introduce the loop model we will use to describe positronium 
decay. 

In that model, the positronium decays into a virtual electron-positron pair 
which subsequently annihilates into real or virtual photons (an odd number for 
ortho-states, an even number for para-states). The coupling of the positronium 
to its constituents is described by a form factor, denoted by Fb- It is not 
assumed to be a constant, since a constant form factor would amount to consider 
positronium as a point-like bound state. Specific forms for Fb will be discussed 
later; for now we just mention that it should somehow be related to the bound 
state wavefunction. 

For parapositronium decay into two photons, our model is represented by 
figure 1. The corresponding amplitude is written 

MT (V-Ps ^F B Tr ^-A—^ J ^ F -} (2) 

with m the electron mass and Fb = Fb (q 2 , P ■ q). The tensor T^ v is the 
scattering amplitude for off-shell e + e _ , with incoming momenta ^P — q and 
\P + q, into two photons : 



T'' 1 ' (e + e~ — > 77) = ie^ = ieY + ieY = ie^ 



The decay width is expressed in terms of this amplitude 

r (jp-Ps -+ 2 7 ) = 1 JL J d ^ J2 \M» V (P-Ps - 2 7 ) e^e* 2u 

pol 

where M < 2m is the positronium mass. 

This model is readily extended to any final states by changing the scattering 
amplitude T to the proper one, and to orthopositronium states by the replace- 
ment 75 -^fi with e the orthopositronium polarization vector. Therefore, the 
conclusions reached in the next section will be of complete generality. 
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3 The Standard Approach as a Dispersion Relation 



What we intend to show in this section is that the formula @ leads to the 
expression for the decay amplitude found in the literature (see for example Q , 

§, I) 

MivPs^H- J , 2 t k 2E J (k 2 )Tr{(l + 7 °) 75 r^ (k,k',h)}e^e* 2u 

(3) 

if we define the form factor as Fb = Cij) (k 2 ) (k 2 + 7 2 ); i\) (k) is the bound 
state wavefunction and 7 2 = m? — M 2 /4 (related to the binding energy through 
Eb = M — 2m ~ — ma 2 /4). The scattering amplitude describes the process 
e _ (k) e + (k') — ► 27 with on-shell electron and positron of momenta k and 
k' = — k and energies E^ = E^/ = Vk 2 + m 2 . 

An interesting feature of this formula is that energy is apparently not con- 
served, since the electron and positron have energies E^ > M/2. The justifica- 
tion of this formula as a simple convolution between the positronium wavefunc- 
tion and the scattering amplitude for e + e~ — > 27 is therefore inadequate. What 
we are going to demonstrate now is that one should understand (3) as a disper- 
sion integral along the loop model branch cut. Further, the appearance of the 
spin wavefunction of the bound state (l + 7 ) 75 is questionable. Indeed, it is 
clear that it is the moving electron-positron pair which has to be projected onto 
the required spin state, so the projector cannot be simply (1 + 7 ) 75. This is 
confirmed by the fact that (||) will be found by neglecting some momentum |k| 
dependence in the exact formula (g). Taking into account the exact projector 
will introduce some forgotten corrections. 

Let us now demonstrate that by using a dispersion relation (see ||, ||) to 
express the loop integration of (^), we will reach (Q). Let us emphasize once 
again that the whole discussion of this section is readily extended to any para- 
or orthopositronium decay amplitude. 

We first compute the imaginary part of (§), ImTy-j (-P 2 ) = Im M.^ u (p-Ps — > 27) £1^2, 
for an arbitrary initial mass P 2 . Considering the two possible cuts (figure 2), 
we obtain Im Tfi by replacing the two propagators on each side of by delta 
functions 

xTV {75 (i - \ + «j) r<" U + ^ + m\ \ e;„e;„ 

After a straightforward integration over q° and |q|, with P = (y/ P 2 ,Q^j, we 
reach 

Im7 >'( p2 ) - 

xTr | 75 (V - ^ + m) (4 + t + rr^j } e* lfl e* 2l/ 
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In the course of the derivation, the delta functions forced q° = and |q| 
y/p 2 /4 _ m 2 j n other words, the electron momenta are 

1 _ I i . . / j „ \ ' , ,■■> 2 




-I'l r/ ~ | \ — • <| I willi (-P±r/) -_-|q|-- m - (I) 

This kinematics is to be understood in the trace evaluation. The angular de- 
pendence arises from the relative orientations of q and the photon momenta 
l\. Note also that the relation (4) cannot be satisfied for the physical value 
P 2 = M 2 < 4m 2 . This is obvious since the loop cannot have an imaginary 
part for the physical bound states, its constituents being always off-shell. From 
the kinematics (Q) one can prove that the factors on both sides of T^ u are true 
projectors, which serves to enforce gauge invariance in the expression 

^ + m) ^[i+^P + m 

Indeed, those two projectors play exactly the same role as external spinors when 
demonstrating Ward identities. Remark also that the dependence of the form 
factor on the loop energy can be taken arbitrarily since q° = 0. 

The real part will now be calculated using an unsubstracted dispersion re- 
lation 

where it is understood that P 2 should be replaced by s everywhere, i.e. scalar 
products that will appear when evaluating the trace should be expressed with 
the kinematics defined for an initial energy s. Since M 2 < 4m 2 , the principal 
part can be omitted and 7jj (M 2 ) = Re7/j (M 2 ). Now let us write the form 
factor in the general form 

F B = C^T (q 2 ) (q 2 + 7 2 ) = C<j) T (a/4 - m 2 ) • (s - M 2 ) /4 (6) 

with 7 2 = m 2 — M 2 /4 and 4> the bound state wavefunction at zero separation. 
Then Eq. (5) can be written as 

m 1 



64vr 2 



x Tr | 75 [i - £ + m) T^ (i + £ + m) } e\^ 2v 

Let us transform the s integral back into a |q| integral, keeping in mind the 
constraints obtained when extracting the imaginary part. Using q 2 = s/4 — m 2 , 
ds = 8 |q| d|q|, the decay amplitude dispersion integral is 



where, as the notation suggests, it is understood that any P 2 appearing in 
the amplitude must be replaced by 4|q| 2 + 4m 2 . In particular, \JP 2 (q) can 
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be replaced by 2£q with = y|q| 2 + m 2 . This amounts to consider the 
scattering amplitude with incoming on-shell electron-positron having momenta 

(j}P (q) ± qj = m 2 (since q° = 0). Note the fact that £q > M/2, apparently 
the energy is not conserved. This is not surprising since the present formula is 
a dispersion integral, done along the cut where P 2 (q) > 4m 2 . Finally, in view 
of the kinematics, we introduce k = \P (q) + q and k! = \P (q) — q (hence 
Ek = Eft = Eq and k = — k' = q) to write the amplitude simply as 



Tji (M 2 ) =^J ff^ [^T (k 2 )] Tr { 75 (- ft + m) I*" (fc, fe', fc) + m)} e^e 

(7) 

where r^" (fc, Zi) is the amplitude for on-shell e (k) e + (k') scattering into 2j. 
Gauge invariance is present due to the two projectors, well defined since k 2 = 
k' 2 = m 2 . Hereafter, the equation (7) will be referred as the standard approach 
expression. Indeed, we can recognize (|7|) as the standard decay amplitude for 
bound states || : 

, r f]3lr 1 1 

M {p-Ps 27) = V2M J ^3 J, (k) -j==-j==M (e~ (k) , e + (-k) - 2 7 ) 

(8) 

with the amplitude constrained to the required spin state given by 

M (e- (k) , e+ (-k) -» 2 7 ) = ^-|^rr {75 (- ^ + m) (fe, fe', Qi + m)} ej^ 

Matching the present expression for the amplitude (7) with the expression (||), 
we obtain the constant C 

C = y/M/m (9) 

and the form factor T is identified with the wavefunction as (p J- (k 2 ) = ip (k). 
This means that the function T (k 2 ) is normalized to unity and behaves as a 
delta of the momentum in the limit of vanishing binding energy : 



1v 



/|H k2) = 1, S^( k2 H 2 ' )3i<3,(k) < 10 > 



From the expression (7), it is clear that (|3|) is an approximation. Indeed, 
neglecting k dependences in the two projectors (—fi' + m) and Q6-\-m), we 
reach (|3|) from (R). The projectors appearing in ([?]) for particles in motion will 
introduce some new corrections to the positronium decay rate of the order of 
the binding energy 7 2 , i.e. a 2 . Further, when generalizing to orthopositronium 
decay amplitudes, the expression (0) is gauge invariant while (||) is not. The 
fact that for the parapositronium decay into two photons @ is gauge invariant 
is due to the peculiar feature of 7 s appearance in the trace (see (pi]) below). 

In conclusion, the equivalence between the loop model and the standard 
expression (|?j) for positronium decay amplitudes opens new possibilities for 
explicit computations. This is precisely what we are going to exploit in the 
following sections. 
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4 Form Factor Dispersion Relation 



We have established the correspondence between (g) and (^). Let us construct 
an alternative, but equivalent, dispersion procedure specific to the two-photon 
case that will be used in explicit calculations. By computing the trace in (||), 
the tensor structure factorizes 

M (p-Ps -> 77) = 8me 2 e^hJ 2 ^el^l (m 2 ) (11) 

In this equation, X (M 2 ) can be viewed as an effective form factor, modeled 
as the electron-positron loop with the coupling F&. There is only one term in 
X (M 2 ) since the direct and crossed amplitudes are equal under q — > —q, i.e. an 
allowed variable change as Fb (q 2 , P • q) = Fb (q 2 , —P • q) , and we write 



A r, 1 1 1 

,2 



( 27r ) (l-k P ) -m 2 (q + \p) -m 2 (q-\P + h 

(12) 

It is to evaluate the effective form factor X (P 2 ) that we will now use dispersion 
techniques. Remark that the factorization of the tensor part is interesting, since 
gauge invariance is manifest, and that X (P 2 ) is convergent while the amplitude 
(2) is superficially divergent. 

The factor rj is introduced because there is a subtlety in the above factoriza- 
tion. Indeed, there is an arbitrariness in the choice of variable for the dispersion 
integral. This situation is well-known for the photon vacuum polarization : 

(fc 2 ) = (k 2 g^ - Wk v ) IT (fc 2 ) = (V - ^) n 2 (fc 2 ) (13) 

where one writes a dispersion relation for IT (fc 2 ), which is less divergent than 
II2 (fc 2 ) due to the factorization of the tensor structure. The analogue of (13) 
here is 

M (p-Ps -» 77) = Sme 2 ^?!,,!^^^ (m 2 ) = Sme^^j^el^X, 

Here we will choose to write a dispersion relation for X2 (M 2 ), and this corre- 
sponds to the choice r\ = P 2 /M 2 . This choice seems arbitrary, but is in fact 
necessary to recover the results of the first section, i.e. the standard expression 
(H). To understand this, consider the factorized amplitude (|il~l). In the disper- 
sion procedure used to recover the standard approach in the previous section, 
the dispersion relation was built on the whole amplitude, so it is clear that the 
photon momenta appearing in the tensor structure were also incorporated, i.e. 
they were reduced as l\ )P —>■ l\ tP x VP 2 /2. That is the reason why we must 
include a factor P 2 into the effective form factor X (P 2 ) ■ 

Let us give a general expression for X (P 2 ) as a dispersion integral. The 
loop integral X (P 2 ) has an imaginary part obtained by cutting the propagators 



M 



2 N 



lmX(P 2 " ' P 



1 , d A q 2m5((q-\p) 2 -m 2 ^j27Ti5^(q+lP) 2 -m 
2 M 2 J J2^ FB (q-lp + iy - 2 
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F B 



y/l - 4m 2 Jf 



8vrM 2 J 4vr 1 + yT - 4m 2 /scosi 



1 (s - 4m 2 ) 



(14) 



with s = P 2 . The form factor is evaluated for |q| 2 = s/4 — m 2 and g° = 0. The 
unsubstracted dispersion relation is 



Re J ( M' 



P r +ao ds 



TT 



M 2 



Iml (s) 



(15) 



In our domain M 2 < Am 2 and 1 (M 2 ) = Rel(M 2 ). With the change of 
variable |q| 2 = s/4 — m 2 , ds = 8|q| d|q|, we write 



l(M 2 ) = Rel(M 2 



C<f> f d 3 



'q 



2M 2 J (2k) 



T q 



q| 2 + m 2 + |q| cos i 



(16) 



where F B = C<\> T (q 2 ) (q 2 + 7 2 ) (see 6). 

The final expression of the amplitude is therefore 



M(p-Ps^ 11 )=8me 2 ^ 

and the decay rate is expressed as 

r / T3 \ 47ra 2 m 2 ^ 2 , 1 |2 

r (p-Ps -» 77) = — — 



Aivpu) 7 * * 



rf 3 



q 



(27T)- 



q 



1 



q| + m 2 + |q| cos i 
(17) 



(i 3 q 



■t q 



q| 2 + m 2 + |q| cos i 



(18) 

This is our third representation for the same decay amplitude : the first is 
the loop integral (2), the second is the well-known amplitude (0) or (JsJ) (with no 
approximation for projectors) viewed as a dispersion integral for the amplitude, 
and the third is the present dispersion integral (|l7]) for the effective loop form 
factor X (P 2 ) ■ All three procedures are strictly equivalent to each other. 



5 Positronium Decay to two Photons 

Before going through the rate calculation, we will analyze the delta limit for 
the form factor as in (10). Then we will go through two different calculations 
of r [p-Ps — ► 77), obtained for specific choices of Fg (or equivalently, T (q 2 ))- 

5.1 Decay Rate in the Static Limit 

To compute the decay rate in the limit 7 2 — > for the form factor, we do not 
need to specify Fb- We just need to know that 

F B = C4>oF (q 2 ) (q 2 + 7 2 ) 7 -° + C^ (2VT) 3 5^ (q) (q 2 + 7 2 ) 



s 



By setting T (q 2 ) = (2tt) 3 6^ (q) in (12), we must recover exactly the lowest 
order decay rate T (p-Ps — ► 77) = ^a 5 m. Using (18), we get 



d 3 c 

q| + m 2 + |q| cos 9 



m 



Importantly, this result is independent of the binding energy : the loop do not 
introduce any corrections in the static limit. The decay rate in that limit is 
therefore : 

T(p-P S ^ 11 ) = la 5 m(^C 2 ) (19) 



with \4>o\ 2 = a 3 m 3 /87r. It remains to match C such that purely kinematic 
corrections vanish (M factors in the above formula arise from products like 
(Zi • I2) = M 2 /2 and from the 1/2M decay width factor, while m comes from 
electron propagators in the loop and from the wavefunction (fto). With the 



definition (9) C = \JM /m, the decay rate is exactly V (p-Dm — > 77) = m 
as it should. In other words, the value for C obtained by matching (0) and (Q) 
is such that no correction arise from factor M/m in the static limit. 

To conclude, let us repeat that we have not specified the form factor. This 
means that any form factor which has a three-dimensional delta function limit 
for 7 2 — > gives the correct lowest order decay rate \oC 3 m. In the following, we 
shall present two forms, both built on the Schrodinger momentum wavefunction. 

5.2 Schrodinger Form Factor 

We can now apply the formulas of the preceding section to write down the 
dispersion integral for the form factor 



8 ' 7 (20) 



(q 2 +7 



with 7 2 related to the binding energy and the fine structure constant through 
Eb = M — 2m = —ma 2 /4. This form factor is just the Shrodinger momentum 
wavefunction for the bound state (note that it satisfies the properties (10) as 
used in j|]). 

By using the formula (14) and (6), the imaginary part is after the angular 
integration 



Iml(s) = n x 

V ' 16vrM 2 



87T7 



(q 2 +7 2)2 



2 1 2 

q +7 



xln 

- q 2 =s/4— m 2 



1 + II _ 4m£ 



1 _ /I _ 4m£ 



We now integrate the imaginary part using the unsubstracted dispersion rela- 
tion (|l5|). As always, M 2 < Am 2 so that the principal part can be forgotten 
and I (M 2 ) = ReX (M 2 ) . For the given form factor Ti (q 2 ) , the calculation of 
integral X (M 2 ) is now straightforward, and we get 



J (m 2 ) 



M 3 7T 27 ^ ' 
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The integral needed for this calculation is 



dx 
o x Q - x 



In 



1-VT 



= 2 arctan 



V x o ~ 1 



(22) 



and its derivatives (d /dx ) n . With the result for I, the decay rate is 



2 MY 
— arctan — 

7T 27/ 



1 / 4?77. 2 

T (p-Ps 77) = -a m I 



where we used \4>q\ = a 3 m 3 /87r and C = \[Mjm. By expanding this re- 
sult around 7 = 0, and expressing corrections as a series in the fine structure 
constant a, we recover the standard result as zeroth order : 



V ^ UJ 2 V vr 8 96vr 64 V 

w ^a 5 m (l - 0.637a + 0.351q 2 - 0.166a 3 + 0.074a 4 + O (a 5 )) 



Numerically, the corrections at the one-loop level are 



^Form Fact. 



(p-Ps -» 77) = r o x 0.9954 ps 7.9956 x 10 9 sec -1 



If we combine the present correction with radiative corrections up to order 
a 2 In a as found in the literature [2] , the theoretical value is modified as 

T Rad.Corr. + Form Fact. ^_ pg ^ ^ = ? x ^9^-1 

while the experimental value is (|l|) . We can therefore note that the present form 
factor leads to a too small value for T (p-Ps — > 77), since it introduces new order 
a corrections. The problem is in the form factor, which do not converge fast 
enough towards the static limit delta. In other words, for a given 7 2 , Ti (q 2 ) 
is not enough peaked around q = 0. 

5.3 Improved Schrodinger form factor 

The second possibility we analyze is 

327T7 3 



(q 2 + 7 2 ) 3 



(23) 



This form factor has the same delta limit property (10), but converges faster 
than (|2"C)|). We can therefore expect much smaller corrections than with Ti (q 2 ). 
Repeating the derivation of the preceding section, we get from (113) : 



Iml is) 



C(f>o 

16vrM 2 



32tt7 3 



(q 2 + 7 2 ) 



2\3 



2 1 2 

q +7 



xln 



- q 2 =s/4— m 2 



4m 2 



4m 2 
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The unsubstr acted dispersion integral is (|l5|), and we get using the result (B^) 



j(M 2 ) = ReT(M 2 ) 




M 

arctan — 
27, 



The rate is then 



T (p-Ps — > 77) 



5 

-a m 



4m 
M 2 " 




arctan • 



M 
27, 



Thus, transcribing into a series in a : 

1 /" 1 2 T 
T (p-Ps -> 77) = -a 5 m ( 1 - -a 2 H a 3 -a 4 + O (a 5 

2 \ 4 3-7T 64 V 

Numerically, the corrections starting at order a 2 are 

r (p-Ps -> 77) = (1 - 1-32 x icr 5 ) x r„ 

i.e. very small. If we combine the present correction with radiative corrections 
up to order a 2 In a as found in the literature, the theoretical value is modified 

as 

r rad.Corr. + Form Fact. ( p _ pg ^ ^ = 7>ggg4 x 1Q 9 ^-1 

As announced, this form factor leads to acceptable corrections, contrary to the 
J- 1 (q 2 ) form factor. 

Before closing this section, let us discuss how the form factor Tu (q 2 ) could 
arise as the right form factor. First, one can see that it implies the Fb form 



F B = Ccf> c 



32tt7 3 



(q 2 +7 



2^3 



(q 2 + 7 2 ) =C4> 



(q 2 + 7 2 ) 2 



4Y 



The factor in brackets is the momentum Schrodinger wavefunction while the 
additional factor 7 2 gives to the coupling Fb the desirable property of vanishing 
when 7^0. From known Bethe-Salpeter analyses (see for example 
one constructs an approximated bound state wavefunction by considering only 
coulombic photon exchanges in the Bethe-Salpeter kernel. Let $ (q) be this 
wavefunction (the Barbieri-Remiddi wavefunction J?J), but with its spin wave- 
function part omitted. Then, if one defines the coupling as 



F B = C x 4 7 2 x (q 2 + 7 2 ) x * (q) 



(24) 



by going through the dispersion analyses of section 3, one ends up with ip (q 2 ) = 
J- 1 1 (q 2 )- The reason for this is quite technical, but let us just mention that 
usually one "uses" some part of \P (q) to set the energy q° to zero to go from 
a four-dimensional towards a three-dimensional convolution- type decay ampli- 
tude (see [3], [£|), while here it is the dispersion relation which enforces q° = 0, 
leaving \& (g) unaltered. This brings a supplementary 1/ (q 2 + 7 2 ) factor. Fur- 
ther, the formula (|24| ) is the analogue of standard expressions making the con- 
nection between Bethe-Salpeter vertex and bound state wavefunction via some 
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propagators. Even if the present remarks do not constitute a rigorous proof, 
they point towards Tn (q 2 ) as the appropriated form factor, rather than the 
widely used J-j (q 2 ). In conclusion, we would like to stress that the well-known 
Bar bier i-Remiddi wavefunction, used as a basis for decay rate calculations, is 
not simply the Schrodinger momentum wavefunction, contrary to what is often 
seen in the literature. 



6 Conclusions 

The result of our analyses is three-fold. First, we have shown how the standard 
convolution- type factorized amplitude can be given a coherent grounding from 
dispersion relations. This demonstrates that the formula usually quoted is an 
approximation, missing some of those O (a 2 ) corrections it is meant to evaluate. 
Then we have calculated those corrections for p-Ps — * 77 using our exact 
factorized formula, or equivalently, its effective form factor realization, for two 
different couplings Fb- Finally, we comment on those couplings. 

Concerning this last point, it should now be clear that the dispersion method 
used to factorize the bound state dynamics from the annihilation process im- 
poses some constraints. By this we mean that there is a definite procedure for 
gluing together the bound state wavefunction and the decay amplitude. 

More precisely, the energy dependences and the spinor part of the Bethe- 
Salpeter wavefunction are usually incorrectly introduced. As we have discussed, 
the energy dependences introduce a modification of the naive Schrodinger wave- 
function, while the spinor part should be introduced in the decay process rather 
that in the bound state wavefunction, in order to properly project the con- 
stituents onto the required spin state. 

Finally, we have shown that the static limit reproduces the first approxi- 
mation to decay rates. This approximation will be modified by higher order 
corrections only, provided the form factor has a three-dimensional delta limit 
for 7^0. Whether this limit property is shared by quarkonia wavefunctions 
is still an open question. 



In a following paper, the present formalism will be applied to the orthopositro- 
nium decay into 777. The analytical defects of this decay channel, noticed in 
|1C], will be completely removed thanks to dispersion techniques. As a result, 
in addition to the new O (a 2 ) corrections found here, a whole set of new am- 
plitudes contributing to that order will be presented. Those amplitudes are 
non-factorizable contributions to the decay rate that could not appear in the 
standard factorized procedures, but unavoidable from the requirements of gauge 
invariance and analyticity. 
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Parapositronium Decay to yy : Figures 




Figure 2 : The two cuts contributing to the imaginary part 



